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ABSTRACT 

We conduct a statistical analysis of the radio source population in galaxy clusters as a func- 
tion of redshift by matching radio sources from the Faint Images of the Radio Sky at Twenty- 
Centimeters (FIRST) catalog with 618 optically-selected galaxy clusters from the first Red- 
Sequence Cluster Survey (RCS1). The number of excess radio sources (above the background 
level) per cluster is 0.14 ± 0.02 for clusters with 0.35 < z < 0.65 and is 0.10 ± 0.02 for clusters 
with 0.65 < z < 0.95. The richest clusters in the sample have more radio sources than clusters 
with low or intermediate richness. When we divide our sample into bins according to cluster 
richness, we do not observe any significant difference (> 1.5c) in the number of radio sources 
per unit of cluster mass for the galaxy clusters with 0.35 < z < 0.65 as compared to the galaxy 
clusters with 0.65 < z < 0.95. Thus the entire sample can be characterized by the number of 
(Liaghz > 4.1 x 10 24 W Hz -1 ) radio sources per unit (10 14 M ) mass, which we measure to 
be 0.031±0.004. We further characterize the population of galaxy cluster-related radio sources 
through visual inspection of the RCS1 images, finding that although the radio activity of bright- 
est cluster galaxies (BCGs) also does not strongly evolve between our high and low redshift 
samples, the lower-redshift, richest clusters are more likely to host radio-loud BCGs than the 
higher- redshift, richest clusters or poorer clusters at the 2-ct level. 

Subject headings: galaxies: clusters: general — galaxies: active 



1. Introduction 

Historically, observing the environments of 
powerful radio sources has been a useful method 
of finding high- redshift galaxy clusters (e.g., 
Dickinson! |l994j). Also, much detailed study 
has been done of the radio sou rce population 
in nearby galaxy clusters (e.g., Ledlowl Il99 

enll 



Owen. Ledlow fe KeeJl995HLedlow fc Owenll99 
Now that uniformly selected samples of more dis 
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tant galaxy clusters exist, we can address the red- 
shift evolution of the population of radio sources 
in galaxy clusters; that is, given a complete sample 
of galaxy clusters, how many of them contain pow- 
erful radio sources as a function of cosmic time? 
The answer has implications for galaxy formation, 
galaxy cluster physics and the selection function 
of large SZ-selected galaxy cluster surveys. 

Heating from powerful radio jets of AGN is 
one feedback mechanism invoked to explain the 
observed X-ray properties of the centers of galaxy 
clusters (Bin nev fe Taborll995[|McNamara fe Nulsen 



20071 f or a r eview) Studies suc h aslBurnsTiji 990f ). 
Eilekl (|2004f ) and iMittal et ail (|2009h have pro- 



vided evidence that most (and maybe all) galaxy 
clusters with X-ray cooling cores also contain pow- 
erful radio-loud AGN in their brightest cluster 
galaxies (BCGs), suggesting that the radio ac- 
tivity of the black hole of the central galaxy is 
linked to the thermodynamics of the cluster gas 
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near the core. Thus, understanding how the pop- 
ulation of radio-loud AGN is related to cluster 
mass and how this population evolves with red- 
shift is important to understanding the energy 
budgets in the gaseous cluster cores. Radio-loud 
AGN feedback may also be important in galaxy 
evolution. For example, applying semi-analytic 
mod els to the output from the Millenium Simula- 
tion, Croton et al. ( 20061 ) found that not only do 
radio AGN solve the problem of the over-cooling 
of gas in the cores of galaxy clusters, but that the 
feedback from radio AGN affects the colors, mor- 
phologies and stellar ages of the massive galaxies 
in clusters, bringing them into better agreement 
with observations than models without such feed- 
back. 

Current experiments such as the South Pole 
Telescope and the Atacama Cosmology Telescope 
use the Sunyaev Zel'd ovich (SZ) effect to dis- 
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galaxy clusters ( Vanderlinde et al. 201dt 



Menanteau et al. 2010l) . These experiments oper- 
ate at millimeter and sub-millimeter wavelengths 
and take advantage of the SZ effect's nominal in- 
dependence from cluster redshift. While most 1.4 
GHz radio sources have little power at such high 
frequencies, a population of radio sources with flat 
or inverte d spectra could contaminate these sur- 
veys (see iLin et al.ll2009l for an investigation of 
the spectral properties of cluster radio sources). 
Little is currently known about the redshift evo- 
lution of the radio sources associated with galaxy 
clusters, but it could affect the cluster sample se- 
lection as a functio n of redshift for these surveys 
(|Lin fe Mohrll2007l ): investigating the redshift evo- 
lution of the radio source population in clusters 
even at 1.4 GHz provides a useful constraint to 
high frequency experiments. 

In this paper, we study the population of ra- 
dio sources from the Faint Images of the Ra- 
dio Sky at Twenty-Centimeters (FIRST) survey 
that lie within the projected inner 0.5 Mpc of 
galaxy clusters fro m the first Red-Sequenc e Clus- 
ter Survey (RCS1; iGladders fc Yeell2005l) . The 
optically selected RCS1 galaxy cluster sample is 
well-characte rized over a wide ra nge in redshift 
(0.35 to 0.95: IGladders et al.lE)07h . allowing us to 
systematically address the evolution of the galaxy 
clusters' radio source population using a uniformly 
selected sample. 

Throughout this work we assume a flat ACDM 



2. Data 

2.1. Red-Sequence Cluster Survey 

The RCS1 data cover a total of about 90 square 
degrees in the Rc and z' filters. The data were 
taken at the Canada-France-Hawaii Telescope 
(CFHT) and the Cerro Tololo Inter- American Ob- 
servatory (CTIO) 4m. IGladders fc Yeel (|2005h de- 
scribe the details of the data acquisition, data re- 
duction and galaxy cluster survey generation and 
characteristics. The survey is designed to take 
advantage of the red sequence cluster detection 
algorithm (|Gladders fc Yeell2000l ). 

In this sample, galaxy cluster centers are deter- 
mined by the center of light from red galaxies, 
and redshifts for the galaxy clusters are com- 
puted using observed red-sequence colors. The 
redshifts are typically accurate to better than 
0.05 for 0.2 < z < 1.0, as indic ated by follow- 
up spectroscopy (B i ndert 20061) and m odeling 



of the RCS data ( Gladders et all l2007h . See 



Gladders fc Yeel (|2005j) for details on the deriva- 
tion and uncertainties on these properties. The 
cluster richness parameter B ffc is used as a mass 
proxy and is defined as the amplitude of the 
galaxy— cluster center correlation function, which 
scales as the number of galaxies within 0.5 Mpc 
(evaluated for h = 0.5) of the cluster center nor- 
malized by the gal axy luminosity functio n and 
spatial distribution. lYee fc Ellingsonl ( 2003 ) show 



that B 9C scales well with other mass proxies using 
the CNOC cluster sample. They empirically de- 
rived the following relation based on dynamically 
determined mass measurements in order to con- 
vert from richness to cluster mass, which we have 
scaled to give us M 2 oo for h = 0.705: 

logM 200!h=0 . 705 = (1.64±0.28)logS gc +(9.90±0.89) 

(1) 

We use B SCi fl (calculated with h = 0.5), which 
is derived using galaxies in the red sequence. B gCv R 
is expected to trace the cluster mass over a wider 
redshift range than B sc , as the red galaxies evolve 
slowly and the cluster red sequence was well in 
place by z = 1. ^> gc ,R is less likely to be affected 
by variations and evolution in the blue fraction 
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of cluster galaxies. Indeed the cosmological anal- 
ysis of the RCS1 survey does not detect evolu- 
tion in the mass-rich ness relation (althoug h the 



constraints are weak; iGladders et al.ll2007i ). Be- 



cause the above mass-richness relation was derived 
for B gc rather than B sc #, we have applied a cor- 
rection factor of 1.23 to the TS gc ,R measurements 
when we calculate mass from richness. We empir- 
ically determined this factor by using both B SCv r 
and B gc for the 212 clusters in our sample which 
overlap with the redshift range of the CNOC clus- 
ters (0.35 < z < 0.54) and ca lculating the mean 
ratio o f the two richnesses. See Yee fc Lopez-Cruzl 
(|l999h or a comparison between Abell richness 



class and B gc (the median B gc ,R of a cluster with 
Abell richness class is about 1090 for h — 0.705). 

We impose a richness cut so that only galaxy 
clusters with B gc R > 300 are included. We also 
require that the uncertainty on the richness be less 
than half the richness value and that the signifi- 
cance of each cluster detection be greater than 3.3 
a. This significance limit was applied in the cos - 
mological study of the RCS (|Gladders et al.ll2007l ). 
which found that for clusters above this limit, the 
contamination rate is < 0.5%. The redshift range 
for the clusters is restricted to 0.35 to 0.95. The re- 
sulting subsample of galaxy clusters is well-defined 
within a cosmological context over this broad red- 
shift interval, as we are using the same sample 
constraints us ed in the cosmological analysis of the 
RCS1 survey (|Gladders et al.ll2007l ). 

The combined area of overlap between the 
RCS1 and the FIRST surveys is about 40 square 
degrees. We limit the cluster sample so that it 
contains only galaxy clusters that are within and 
further than 38' from the edge of the FIRST sur- 
vey, which is just greater than 10 Mpc at a redshift 
of 0.3. Because we also limit our cluster sample to 
clusters with z > 0.35, the 10 Mpc area surround- 
ing each cluster in our sample fully lies within the 
footprint of the FIRST survey (thus, we may use a 
10 Mpc annulus around the clusters' centers when 
calculating the background level of radio sources, 
see section 13". 1.1[) . 

Constraining the RCS1 cluster catalog in this 
way yields a sample of 618 galaxy clusters. 



2.2. FIRST Radio Survey 

The FIRST survey was conducted at the Na- 
tional Radio Astronomy Observatory (NRAO) 
Very L arge Array (VLA) bv lBecker. White. &: Helfand 
(jl995l ). The survey consists of observations at 1.4 
GHz. The VLA was in its B configuration, giving 
the resulting maps a resolution of 5". The survey 
reached a point source sensitivity depth of ~1.0 
mJy. Radio s ource s were selected and cataloged by 



White et al.l (|1997l ) . Source positions are accurate 
to better than 1", enabling matches to individ- 
ual optical objects in many cases. Because of the 
high resolution of the survey data, radio sources 
that are extended, particularly radio sources with 
double lobe morphologies, are often resolved and 
counted as multiple radio sources. We address 
these multiple component radio sources in section 
15X41 

3. Analysis and Results 

3.1. Statistical Analysis of Radio Sources 
around Clusters 

3.1.1. Methods 

Radio luminosity is calculated from flux using 
the following relation: Li AGHz — ^Fiaghz Da 2 (1+ 
z) 3+Q , where we assume a spectral index of 
a = 0.8 where L u oc v~ a . According to this rela- 
tion, a radio source at the flux limit of the FIRST 
survey (1 mJy) and at the redshift limit we are im- 
posing on the RCS1 survey (z = 0.95) would have 
a luminosity of 4.1 x 10 24 W Hz -1 , so we adopt 
this as our luminosity cutoff. At this luminosity 
it is likely that radio emission is dominated by 
AGN activity rather than star formation, as radio 
sources dominated by star formation most ly lie be- 
low 10 23 W Hz" 1 dMauch fc Sadler! 12007?) . While 
there is some evidence that this crossover limit 
in the radio luminosity function evolves with red- 
shift, our sa mple is still well above the limit eve n 
at z = 1 ( Owen. Strazzullo. fc Pannella 2009). 
Also, it has been shown in nearby clusters that 
the crossover between the p opulations is fainter i n 
clusters than in the field (jMiller fc Owenl l2002t ) . 
further ensuring that cluster radio sources brighter 
than the luminosity limit we apply are dominated 
by AGN emission. 

For each galaxy cluster in our sample, we iden- 
tify all nearby FIRST radio sources. The clus- 
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ter redshift is used to determine the corresponding 
position (in proper megaparsecs) with respect to 
the cluster center and luminosity for all of these 
nearby radio sources. The luminosity cutoff de- 
termined above is then imposed so that only radio 
sources above this luminosity limit are counted in 
this statistical analysis. We repeat this procedure 
for every galaxy cluster in the sample, compiling 
the list of positions and luminosities for all radio 
sources around all clusters. Adopting the clus- 
ter redshifts for background and foreground radio 
sources does not bias our results in any way be- 
cause we use the resulting positions and luminosi- 
ties of radio sources with respect to the clusters 
to statistically determine the background level for 
the sample. 

Because of the high resolution of the FIRST 
data, extended sources (for example, two lobes 
of an AGN) are split into multiple sources in the 
catalog. These sources are thus counted twice in 
this statistical analysis. Also, if a radio source 
is nearby two galaxy clusters, it will be counted 
twice, once for each cluster. We addressed these is- 
sues by visually inspecting the galaxy clusters that 
contain radio sources within a projected 0.5 Mpc; 
see section [XU While there were 16 galaxy clus- 
ters with centers within 0.5 Mpc of another galaxy 
clusters as determined at the redshift of one cluster 
(7 of which were pairs, 2 of which were within 0.5 
Mpc of other clusters according to their redshifts 
but not according to the other clusters redshifts), 
no radio sources overlapped between two different 
clusters. 

In order to calculate the background level of 
radio sources, we determine the number of radio 
sources above the luminosity cutoff (with the ra- 
dio source luminosities determined for the clus- 
ter redshifts as discussed above) in a 5 to 10 Mpc 
annulus around the cluster centers. We extrap- 
olate this constant background surface density of 
radio sources to the inner cluster regions, subtract- 
ing the background contribution to calculate the 
number of excess radio sources. The uncertainties 
are calculated assuming Poisson statistics for both 
the background and the signal. Figured] shows the 
signal-to-noise ratio in radial bins around the clus- 
ters' centers of the number of excess radio sources 
above the background level. 

As can be seen in the figure, most of the excess 
radio sources are within 0.5 Mpc of the cluster 



center, so we define our statistical sample of radio 
sources in galaxy clusters as the excess number 
of radio sources within a projected radius of 0.5 
Mpc of the clusters' centers. The total number of 
radio sources (including background sources) that 
we find within 0.5 Mpc of clusters' centers is 111, 
which is significantly higher than the background, 
from which we would expect 36.2±0.3. The cluster 
with the highest number of radio sources within 
0.5 Mpc has 4, although 3 of these are multiple 
components of the same radio galaxy. For a dis- 
cussion of the frequency of multiple radio sources 
(both independent and multiple components) in 
clusters, see section 13.2.41 

Correctly handling the background population 
of radio sources is important in this analysis be- 
cause the lower redshift objects subtend more 
area of sky than high redshift objects, so poor 
background estimation could mock redshift evolu- 
tion. We check our background estimation against 
a bootstrap analysis in which we calculated the 
number of radio sources around random positions 
within the area of overlap between RCS1 and 
FIRST. Restricting this analysis to the survey re- 
gion ensures that we did not sample more or less 
sensitive regions of the FIRST survey than our 
data cover so that unevenness in radio coverage 
does not affect our background estimation. We 
assigned each random position a "redshift" drawn 
from the actual cluster redshifts of the galaxy clus- 
ters. The background-subtracted number of ra- 
dio sources per random position is consistent with 
zero. This remains true for both the random po- 
sitions that were assigned redshifts less than 0.65 
and for random positions that were assigned red- 
shifts greater than 0.65. Thus, we see no evidence 
for large scale structures correlated on scales of 
5 to 10 Mpc from the clusters' centers biasing 
the background level we calculate, and we do not 
see any bias from the background calculation that 
could introduce discrepancies between the high 
and low redshift samples. 

3.1.2. Redshift Evolution 

The total number of radio sources within 0.5 
Mpc of clusters' centers in our low redshift (0.35 

< z < 0.65) sample is 62 radio sources in 331 
galaxy clusters, and in our high redshift (0.65 < z 

< 0.95) sample is 49 radio sources in 287 clusters. 
The number of excess radio source components per 
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cluster (after background subtraction) is 0.14 ± 
0.02 for the low redshift sample and 0.10 ± 0.02 
for the high redshift sample. Errors are estimated 
assuming Poisson statistics in the signal and back- 
ground number of radio sources. The significance 
of the difference between the high and low red- 
shift samples is 1.4c. Thus, the level of evolution 
in this sample is constrained to be small (espe- 
cially in comparison with the evolution in cluster 
X-ray AGN, see discussion in section l4"7Tj) . 

An important correlation we notice in our data 
(w hich has be en previously noted, see for exam- 
ple [Best][2007) is that the richest clusters are more 
likely to host a radio source than poorer clusters 
(see Figure [2]). In order to investigate whether 
there are differences in the high and low red- 
shift radio source populations in clusters of similar 
mass, we split our galaxy clusters into three bins 
according to cluster richness (Figured]). The B gc> R 
values of 300, 500, and 800 correspond to masses 
of 1.3xl0 14 Mq, 3.0x10 14 M q , and6.4xl0 14 M©, 
respectively. There are 181 clusters with B SCj # be- 
tween 300 and 500 at low redshift (0.35 to 0.65) 
and 121 at high redshift (0.65 to 0.95), there are 
110 clusters with richnesses between 500 and 800 
at low redshift and 122 at high redshift, and there 
are 40 clusters with B gc ^ 800 and up at low red- 
shift and 44 at high redshift. See Figures |4] and [5] 
for the distributions of the redshifts and richnesses 
of the galaxy clusters. We note that the scatter in 
the mass-richness relation is considerable, so for 
any one cluster the mass is not well measured. 

When we normalize the number of radio sources 
per cluster by the average mass, as derived from 
richness, for each richness bin, the number of ra- 
dio sources per unit of cluster mass is roughly 
constant, with a value of 0.031 ± 0.004 per 10 14 
Mq. Because the masses are derived from optical 
richnesses, this implies that the number of radio 
sources per cluster galaxy is constant, as has bee n 
noted in previous studies ( Ledlow fc Owen 1995). 
Comparing the normalized excess counts of radio 
sources per cluster, we find that the low-redshift 
clusters have more radio sources per unit mass 
than the high-redshift clusters at the 1.5cr level. 
This may be driven by a similar trend in the fre- 
quency of radio activity in the clusters' BCGs (see 
below, section [3.2.2P ; however, the significance of 
the difference is low. Alternatively, the data con- 
strain the redshift evolution to be small. The 1-<t 



uncertainty on the number of radio sources per 
unit mass for each bin in redshift and richness 
ranges from 27 to 53% of the signal, making evo- 
lution by factors of a few unlikely. 

3.2. Visual Inspection 

3.2.1. Methods 

To supplement our statistical analysis, we visu- 
ally inspected images from the RCS1 survey data 
of the galaxy clusters in our sample. We did this 
for a number of reasons: to better investigate the 
relationship between the brightest cluster galaxies 
(BCGs) and the radio sources, to determine the 
frequency of multiple radio sources within single 
galaxy clusters, and to identify radio sources that 
are obviously associated with foreground or back- 
ground galaxies. 

When visually inspecting the cluster images, 
we considered the cluster redshifts and richnesses, 
and we identified bright potential cluster mem- 
bers. Each cluster's center is defined as the center 
of red light, not a BCG position. BCG candidates 
were identified according to the following crite- 
ria: central position in the cluster, morphological 
appearance, position in the color-magnitude dia- 
gram, and the presence of satellites or an extended 
halo of optical light. While in principle we could 
simply select the brightest galaxy within some dis- 
tance from the cluster center and with some color 
near the red sequence, we are more intereseted in 
the central galaxy of the galaxy cluster, and in 
practice the identification of this can be ambigu- 
ous. For example, there may be two bright central 
galaxies where one is slightly brighter, but another 
is more central and has nearby satellite galaxies. 
Or the brightest, dominant galaxy could be offset 
from the cluster center, particularly for clusters 
with extended or disturbed morphologies. Some- 
times multiple candidates were identified as possi- 
bilities. 

Radio sources were also visually inspected. Ra- 
dio sources were matched to optical counterparts 
when possible, and radio morphology was noted. 

As an example of the benefits of this visual 
inspection of the cluster images, Figure [6] shows 
the radio emission in the galaxy cluster RCS1 
J132655+302112. This galaxy cluster is at red- 
shift 0.37 and has a B SCj _r of 1425 (corresponding 
to M 20 o = 1.7 x 10 15 M ). In this example, two 
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components of an extended radio source are identi- 
fied as sources in the FIRST catalog. These radio 
sources appear to be two lobes radiating from a 
bright central galaxy, a good BCG candidate lo- 
cated about 0.1 Mpc from the cluster center (as 
identified by the center of the red galaxies' light.) 
Visual inspection identifies this system, but sim- 
ply matching individual radio sources with indi- 
vidual optical objects using a small radius of tol- 
erance could miss sources like these. This example 
demonstrates the usefulness of characterizing the 
population by eye in order to better understand 
the nature of the population. Of course, both of 
the components of this radio source are within the 
central 0.5 Mpc of the galaxy cluster, so this was 
easily identified in the statistical sample. 

3.2.2. Radio-loud BCGs 

Previous studies have shown that BCGs are 
more likely to be radio -loud than other galaxies 
of similar stellar mass (jBestl 120071 ) , indicating a 
difference in the radio-loud BCG population rel- 
ative to the general cluster radio source popula- 
tion. We investigate whether there is a redshift 
evolution or cluster mass dependence in the radio 
activity of the BCGs in our sample. 

Figure [7] shows the fraction of galaxy clusters, 
divided into bins according to richness and red- 
shift, that contain radio-loud BCGs. Poisson er- 
ror bars are plotted. We find that the fraction 
of radio-loud BCGs does not depend strongly on 
richness for most of the sample. However, the 
prevalence of radio-loud BCGs for clusters in the 
low-redshift, richest bin is higher than the num- 
ber of radio sources per cluster in the other bins, 
possibly indicating evolution in the radio activity 
of the BCGs in the richest clusters. We calculate 
the uncertainty on our measurement of the differ- 
ence between the fraction of radio-loud BCGs in 
the high redshift bin compared to the low redshift 
bin by adding in quadrature the uncertainties for 
each bin. We find that the radio-loud BCG frac- 
tion in the richest low-redshift clusters is larger 
than that in the richest high-redshift clusters with 
a significance at the 1.9<r level. When we com- 
pare the radio-loud BCG fraction for the richest 
low redshift clusters to the radio-loud BCG frac- 
tion of all other clusters, the excess is significant 
at the 2.2a level. The numbers of clusters in the 
richest bins are also comparatively small (9 out of 



40 clusters have radio-loud BCGs at low redshift 
and 3 out of 44 clusters have radio-loud BCGs at 
high redshift.) However, nearby massive clusters 
are the most likely to have short cooling times 
with respect to the Hubble time, potentially al- 
lowing cool gas from the ICM to accrete onto the 
BCG, and possibly onto the central AGN. This 
picture is supported by recent observational stud- 
ies that find a lack of strong cooling core clusters 
at higher redshifts (such as z > 0.5 and z > 0.7, 



Vikhlinin et alj 12007: Sant os et all 12008 ). If this 



scenario holds, then we would expect an excess 
of radio-loud BCGs in the richest nearby clusters 
in comparison with less massive or higher redshift 
clusters, which is what we observe. Further stud- 
ies on larger samples will be necessary to confirm 
or disprove this excess. 

3.2.3. Non-BCG Radio Galaxies 

We identified every radio source whose optical 
counterpart is an early-type galaxy with a color 
consistent with the red sequence of its galaxy clus- 
ter but that is not identified as the cluster's BCG. 
Figure [5] shows the number of such galaxies per 
cluster, divided into bins according to richness and 
redshift, normalized by the mean cluster mass (as 
derived from richness) for each bin. The number 
of non-BCG galaxies in a cluster is proportional 
to its richness, so this is similar to the radio loud 
fraction for these clusters. Some of these clusters 
also contain radio-loud BCGs, but as long as there 
is also an excess of radio-loud galaxies in addi- 
tion to (or instead of) the BCG, the radio sources 
are included in this plot. These cluster member 
candidates were identified based on the color and 
magnitude of their identified optical counterparts 
placing them on the red sequence as well as their 
optical morphology. Unlike in the statistical anal- 
ysis, if two components of an extended source are 
identified as belonging to the same optical coun- 
terpart, the system is counted as one radio source 
in this plot. For each richness bin, the number 
of non-BCG radio sources per cluster (normalized 
by mass) is somewhat higher for the high redshift 
sample than for the low redshift sample. 

Thus, strong evolution in the radio-activity 
in one of those populations (BCG versus cluster 
galaxy) is not masking strong evolution in the 
other population. Also, because we are match- 
ing radio sources with optical sources by eye, we 
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do not count radio sources that are resolved into 
multiple components in FIRST as separate from 
each other. Morphological evolution such that 
the resolution into components is different for the 
high and low redshift samples could mock evolu- 
tion in the radio-loud population in our statistical 
sample. Thus, this visual identification of radio 
sources with BCGs and cluster members upholds 
the overall view that there is little evolution in 
the population of radio sources in clusters from 
redshift 0.95 to redshift 0.35, although the 2-er in- 
dication of an excess of radio-loud BCGs in the 
richest low-redshift clusters is interesting evidence 
for evolution that is worth exploring with a larger 
data set. 

3.2.4. Multiple Radio Sources 

The FIRST catalog splits multiple components 
of extended radio sources into different sources, so 
to explore the frequency of multiple radio active 
galaxies within clusters, it is instructive to view 
them by eye. For example, the radio source in Fig- 
ure E] was identified as two sources in the FIRST 
catalog. We have identified twelve extended radio 
sources in our galaxy cluster sample that are split 
into multiple components by FIRST. These multi- 
ple components contaminate our statistical analy- 
sis, in which we consider radio sources as individ- 
ual and uncorrelated, although the contamination 
is somewhat alleviated because we treat multiple 
components as separate radio sources both in the 
determination of the background level and in the 
determination of the signal. 

Twelve galaxy clusters have multiple radio 
sources projected within 0.5 Mpc from their cen- 
ters that we have identified as distinct from each 
other, two of which have three apparently inde- 
pendent radio sources. Fourteen percent (±2% 
Poisson counting uncertainties) of the galaxy clus- 
ters in our sample have radio sources within 
a projected 0.5 Mpc of their centers (including 
sources we expect from the background level.) Of 
these galaxy clusters, 14±4% have multiple ra- 
dio sources identified as independent sources, and 
17±12% of the sources with multiple radio sources 
identified as independent have three independent 
radio sources. Because all of these percentages 
agree with one another to within their uncertain- 
ties, we note that galaxy clusters with one radio 
source are not much more or less likely to have a 



second radio source, and galaxy clusters with two 
radio sources are not much more or less likely to 
have three. In other words, the radio loudness of 
galaxy cluster members seems to be independent 
of the radio loudness of other galaxy cluster mem- 
bers, although we note that with such a small 
number of clusters in the sample with multiple 
independent sources, the constraints are weak. 

One of the strengths of our statistical analysis 
is that radio sources which are actually lobes of 
a central radio galaxy are included in the sample. 
Because we include all sources within 0.5 Mpc of 
the center of the cluster, we are likely including 
sources of bright radio emission which are associ- 
ated with the BCG even if they are not coincident 
with the BCG because they are lobes. Visual in- 
spection of the radio sources and their relationship 
to the optical galaxies gives us an idea of how im- 
portant this population is. Of the 37 radio-loud 
BCGs we identify in our sample, 9 of them have 
extended radio mor phology. Of these. 6 ar e likely 
FRI's and 3 FRIFs ilFanaroff fe Rilevlll974h based 
on their morphologies. For a recent study com- 
paring the environments of extended and point- 
like radio sources, see (|Wing fe Blantonl l2010h . 
We also calculated the luminosities of these ex- 
tended s ources from their NR AO VLA Sky Survey 
(NVSS; ICondon et al~lll998h fluxes. The NVSS 
was conducted at 1.4 GHz with the VLA in its 
compact D and DnC configurations, so is better 
suited to measuring the fluxes of extended sources 
than FIRST, which m ay resolve out some flux. 
Ledlow fc Owen] (|l996t ) found that for Abell clus- 



ters, the dividing line in the radio luminosity-R- 
band magnitude plane between FRI and FRII ra- 
dio galaxies, which is very sharp in the general 
radio galaxy population (althou gh this sharpnes s 



Lin et all [20101. 



has been recently disputed, see 
is not as distinct in nearby Abell clusters. In our 
sample, all 6 FRI's lie below the line as expected, 
and 2 of the 3 FRII's lie above the line as ex- 
pected, but one FRII is less radio- lum i nous than 
expected. However, Ledlow fc Owen ( 19961 ) de- 
scribe the FRIFs that lie within the FRI part of 
the radio luminosity-magnitude plane as being un- 
usually small, with the lobes confined to the op- 
tical extent of the galaxies that host them. The 
only low-luminosity FRII identified in our sample 
looks fairly typical, with lobes located ~ 200 kpc 
from the host galaxy. Larger samples are needed 
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to better investigate the prevalence of FRII's, both 
luminous and underluminous with respect to the 
FRI/II luminosity division, in clusters at high red- 
shift. 

There were no occasions in this analysis where 
one radio source was within 0.5 Mpc of two cluster 
centers. 

3.2.5. Identification of Foreground and Back- 
ground Sources 

We inspected all of the galaxy clusters that con- 
tain radio sources within 0.5 Mpc of their centers. 
Some of these radio sources are chance projections. 
Given the background estimates from our statis- 
tical analysis, we expect 36.2 ± 0.3 radio sources 
within 0.5 Mpc of clusters' centers by random su- 
perpositions. Radio sources are matched by eye 
with optical counterparts. If the optical coun- 
terparts are elliptical galaxies with colors consis- 
tent with the red sequence at the redshift of the 
potentially associated cluster, the radio source is 
identified as a cluster member. Thirty-two galaxy 
clusters contain radio sources within the projected 
inner 0.5 Mpc that are not identified with cluster 
members. Of these, seven have multiple sources (2 
of the galaxy clusters' multiple sources are iden- 
tified as independent, the others are identified as 
multiple components.) The consistency of these 
numbers provides a check that we are being nei- 
ther too generous nor too miserly in our identi- 
fication of radio loud cluster members, because 
the number of galaxies we leave out as either fore- 
ground, background or indeterminate is consistent 
with what we expect from the background level. 
This consistency also indicates that most of the 
bright radio sources associated with galaxy clus- 
ter members are hosted by early type galaxies that 
lie on the red sequence. 

4. Discussion 

4.1. Comparison with Previous Work 

4.1.1. Radio Galaxy Fraction 

We compare the number of radio sources per 
cluster found in our sample wit h the number of 



radio sources per cluster found in lBranchesi et al 
(|2006l ). who studied radio AGN in X-ray selected 
galaxy clusters with 0.3 < z < 0.8. The FIRST 
survey has a higher sensitivity limit than their 



pointed observations, although they also used the 
VLA operating in the B configuration, so their 
resolution is the same. They also assume a spec- 
tral index of 0.8, and they include radio sources 
within 0.4 Mpc of the clusters' centers. Their clus- 
ter sample is X-ray selected from 87.4 deg 2 of sky 
and contains 18 clusters with 0.3 < z < 0.8, al- 
though only one cluster has z > 0.7. In order to 
compare, we restrict our sample to the most mas- 
sive (B gc ^ > 800) clusters with 0.35 < z < 0.69. 
There are 50 clusters in this richness and mass 
rang e, so we are probing to lower cluster mass 



than lBranchesi et al.l ([20061 ) . who have 18 clusters 



in about twice the area of sky but over a simi- 
lar redshift range, although because their sample 
of clusters is X-ray selected, their sensitivity as a 
function of redshift is different than for the op- 
tically selected sample. When we calculate the 
number of excess radio sources with luminosities 
greater than 1 xlO 25 W Hz -1 within 0.4 Mpc of 
these cluster centers, we find the number per clus- 
ter to be 0.18 ±0.06. According to their Figure 8, 
they find approximately 0.1 ± 0.08 radio sources 
with such bright luminosities per cluster. Thus, 
the studies are in general agreement with each 
other on the number of radio sources in clusters 
in this redshift range. 

In order to compare to an optically-selected 
galaxy cluster s ample, we also compare our data 



to the results of ICroft. de Vries. k Beckerl (|2007l ) 
who correlate brightest cluster galaxies from the 
maxBCG catalog of galaxy clusters in the Sloan 
Digital Sky Survey with FIRST radio sources. 
They use an algorithm to identify radio sources 
with double lobes that are likely associated with 
the BCG, so their radio source sample includes 
both point sources and extended sources. They 
have a much larger sample (over 13,000 clusters) 
but are restricted to low redshift (0.1 < z < 0.3). 
A direct comparison is not straightforward be- 
cause their luminosity cutoff is low (1 mJy at 
z = 0.1 corresponds to about 2.5 x 10 22 W Hz -1 ) 
and the mean cluster N200 richness parameter 



7 , corresponding 


to ar 


(Johnston et al. 


2007 N 



We 

compare the RCS1 clusters that fall between red- 
shift 0.35 and 0.45 with their sample. Using a 
luminosity cutoff of 6.5 x 10 23 W Hz -1 (which 
corresponds to a 1.0 mJy radio source at a red- 
shift of 0.45) we find the number of radio sources 
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per RCS1 cluster in excess of the background 
of radio sources to be 0.29 ± 0.06. We then 
scale this number to account for the differences in 
average mass, in radio luminosities probed, and 
in the fact that we are counting all radio-loud 
galaxies in clusters, not just BCGs. The aver- 
age Bg Ci R of this RCS1 subsample is 627, corre- 
sponding to 4.3 x 10 14 M©. Extrapolating from 
the fiel d 1.4 GHz luminosity fu nction for rad io 



X-ray selected galaxy clusters. They also found no 
evidence for evolution with redshift in the cluster 
radio source population. 

Branchesi et al. ( 20061) calculated the radio lu- 



AGN in iMauch fc Sadlerl (120071) dLedlow fc Owen 
(|l996l ) find no difference between the cluster and 
field radio lumi nosity functions) we expect the 
radio sources in ICroft. de Vries. fc Beckerl (|2007h 
to be about 3.4 times as numerous as the radio 
sources to which we are matching RCS1 clusters 
in this comparison. Finally, from our identifica- 
tion of radio-loud BCGs as discussed in section 
I3.2.2) we find that 75% of clusters with radio 
sources above the level expected from the back- 
ground have radio-loud BCGs. Thus the scaled 
radio-loud BCG fraction for these RCS1 clusters 
is 17.4%. The Poisson uncertainty on this frac- 
tion is about 4%, but the uncertainties in extrap- 
olating across cluster masses and the radio lu- 
minosity function are probably more significant. 
ICroft. de Vries. fc Beckerl (|2007h find that 19.7% 
(±0.4% Poisson uncertainty) of their BCGs are 
radio-loud. Thus without accounting for any pos- 
sible redshift evolution, and with the above ex- 
trapolations for differences in cluster and radio 
source properties, the two studies approximately 
agree. 

4-1-2. Redshift Evolution of Radio Galaxy Frac- 
tion 

Most of the previous work on the evolution of 
radio sources in clusters has been based on in- 
depth studies of relatively small samples of inter- 
mediate to high redshift X-ray selected galaxy 
clusters in comparison to in-depth studies of 
nearby Abell clusters. 

IStocke et al. ( 19991 ) conducted pointed 1.4 GHz 
observations of 19 0.3 < z < 0.8 X-ray selected 
galaxy clusters. They compared the radio activ- 
ity of galaxies in these clusters t o the radio activ- 
i ty ob served in nearby clusters bv lLedlow fc Owen 
( 19961 ) and found no difference between the two 
samples. Exten ding this wor k to p rogressively 
higher redshift, Perlman et al. ( 20041 ) conducted 
pointed 1.4 GHz observations ofl70.5<z<l 



minosity function of 18 X-ray selected galaxy clus- 
ters with 0.3 < z < 0.8. They co mpared thei r re- 
sults to the low-r e dshift studies of Fantil ( 1984h and 
Ledlow fc Oweil (|l996l ) and concluded that their 



radio luminosity function is significantly higher 
than the RLF for nearby galaxy clusters. How- 
ever, as they discussed in their conclusions, their 
X-ray selected clusters have, on average, higher 
mass than the local Abell clusters to which they 
are comparing, and when they take into account 
the presumed di fference in ri c hness between the 
two samples, as Stocke et al.l ( 19991 ) do, they no 



longer discern any difference in the amplitude of 
the radio luminosity function between the high 
and low redshift samples. 

Thus, when the studies of radio AGN in high- 
redshift, X-ray selected clusters properly account 
for the average cluster richness of their high red- 
shift samples, they also see no evidence for evo- 
lution of the radio-loud cluster galaxy population 
with redshift, just as we see no evidence for evolu- 
tion using our uniformly selected sample (see sec- 
tion [3TT2]) . Of course, because these other studies 
use nearby Abell clusters as their low redshift clus- 
ters, they are probing a more local cluster popu- 
lation than we do, and if the radio AGN in clus- 
ters evolve strongly between a redshift of and a 
redshift of 0.3, then these other studies would be 
sensitive to evolution in a redshift range for which 
we have no data. However, the previous studies 
have few clusters with redshifts > 0.6, so we are 
more sensitive to evolution at the higher redshift 
range. 

Best ( 20031 ) conducted deep radio observations 
of a rich z ~ 0.83 galaxy cluster, MS1054-03, find- 
ing 8 spectroscopically confirmed cluster members 
with radio emission brighter than 32/iJy, a large 
excess compared when compared with low redshift 
clusters. Such deep radio observations probe a 
different part of the luminosity function, which 
at low luminosities is dominated by star forming 
galaxies rather than AGN. Instead, we focus on 
the brighter, rarer AGN over a larger statistical 
sample. 
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4.2. Implications 

4-.2.1. The Effect on SZ Cluster Surveys 
generate 



Lin k Mohrl (|2007t) 



predictions of 
point source contamination for SZ surveys. The 
evolution of the radio source population is one of 
their unknowns, so they investigate two models: 
no evolution, and a power-law model in which the 
number of radio AGN in a cluster of mass M scales 
as (1 + z) 7 where 7 = 2.5. For this latter scenario, 
and assuming their propogation of a spectral index 
distribution to determine the predictions of fluxes 
at higher (30, 90, 150 GHz) frequencies from their 
counts at 1.4 GHz, they determine the fraction of 
clusters in an SZ survey that would be lost due 
to contamination from radio sources as a function 
of cluster mass and redshift. Because the SZ flux 
is roughly constant with redshift, and because 
of the increasing luminosity distance with red- 
shift (and with consideration of the fc-correction), 
higher luminosity radio sources are required to 
contaminate the SZ effect from clusters at higher 
redshifts. However, with their power-law redshift 
evolution, they find that for low-mass clusters, 
the evolution of the radio source population wins 
out over these other effects so that the fraction of 
lost clusters is higher at high redshift than at low 
redshift. Conversely, for high mass clusters with 
greater SZ signals, and because of the steepness 
of the radio luminosity function, the fraction of 
lost clusters is higher at low redshift than at high 
redshift. With no evolution, the lost cluster frac- 
tion would be decreased for each redshift bin by 
a factor of (1 + z) 7 , which for the redshifts they 
consider (0.1, 0.6, 1.1) when looking at a range in 
mass (their Figure 15) corresponds to factors of ~ 
1.3, 3.2 and 6.4. 

We do not see evolution in the cluster radio 
source population, and if anything, there is some 
small evidence (1.5er) that lower redshift clusters 
have more radio sources per unit mass than high 
redshift clusters. In order to quantify how signifi- 
cantly our data constrain this type of redshift evo- 
lution, we have performed a chi-square test. We 
calculated the number of radio sources expected 
for (1 + z) 1 evolutionary models, using the clus- 
ter redshifts and normalizing to the total number 
of excess radio sources in clusters that we observe 
in the real data. We compare the models with the 
data for the number of radio sources per cluster for 



the low and high redshift bins (without accounting 
for the richness dependence other than by apply- 
ing a uniform richness cut). The preferred value of 
7 is —1.9 ± 2.8, where the 3tr uncertainties given 
are estimated from fitting a Gaussian to the P- 
values as a function of 7. This test also indicates 
that there is a 0.7% probability that the differ- 
ence between our data and a 7 = 2.5 evolutionary 
model is purely random, implying that such evo- 
lution is unlikely. 

Qualitatively, the main c hange, based on our 



study, in the forecasts from iLin k Mohrl f2007) 
for SZ survey contamination given no evolution in 
the cluster-associated radio source population is 
that the low mass, high redshift clusters are not 
as badly contaminated as low mass, low redshift 
clusters; so, fewer high redshift clusters would be 
lost from SZ surveys. In their more recent study, 
which also includes observations of the spec t ral in- 
dices of radio sources in clusters, ILin et al. (2009) 



adapted a more mild evolution (7 = 1), motivated 
by this work. They found that this makes the 
SZ contamination by AGN much smaller than had 
been previously estimated. 

4-. 2. 2. Comparison with Cluster X-ray and In- 
frared AGN Frequency 

Studies that focus on X-ray selected AGN 
in galaxy clusters observe a strong evolution 
over a redshift range similar t o wha t we probe. 
Martini. Sivakoff. fc Mulchaevl (|2009h observe a 
dramatic increase (by a factor of eight) in the 
population of X-ray selected AGN in galaxy 
clusters over a redsh ift range of 0.05 to 1.3. 
Galametz et al. (2009) conclude that for X-ray- 
selected, optically-selected and radio-selected 
AGN, the enhancement in AGN activity in the 
centers of clusters is greater at 0.5 < z < 1 than at 
z < 0.5. Although the field-corrected radio AGN 
density remains consistent between the redshift 
bins to within their quoted uncertainties, they do 
see some low significance evidence for evolution 
in the radio source population by about a factor 
of two. The clusters in their study are selected 
from optical and infrared imaging data. The evo- 
lution in the X-ray AGN population drives their 
conclusions. 

Our data rule out strong evolution of the radio- 
loud AGN population between a redshift of 0.35 
and 0.95 in our sample. The two AGN populations 



10 



(radio and X-ray) may be evolving differently from 
each other in the cluster environments. Evidence 
of evolution in the X-ray AGN population does not 
impl y evolution i n the cluster radio AGN popula- 
tion. lHart et al.l (|2009t ) identified X-ray AGN and 
radio AGN in 11 galaxy clusters at 0.2 < z < 0.4 
selected from X-ray data, as clusters that are sim- 
ilar to what the Coma cluster would have been 
like at that redshift. They found that X-ray point 
sources and radio sources rarely overlap in their 
sample. The two populations are separate, and so 
predictions that rely on the radio source popula- 
tion evolution (such as projections for the contam- 
ination of SZ surveys) should not use the observed 
cluster X-ray AGN evolution. 

Studies of X-ray AGN in clusters typically de- 
fine a cluster galaxy AGN fraction as their observ- 
able quantity. They typically exclude the BCGs 
from this AGN fraction to avoid pos sible contami 



nation from cluster cool cores (e.g., iMartini et al 



20061 ; I Arnold et ai1l2009l ). Of the different results 
for radio sources in clusters that we have pre- 
sented, the quantity that we measure that com- 
pares most directly to this AGN fraction is the 
number of radio-loud non-BCG cluster members 
per unit mass (as shown in Figure [5]). While 
the slight increase at high redshifts is not signif- 
icant (the significance of the difference between 
the mean number of non-BCG radio sources per 
cluster per unit mass for the high redshift bins 
and for the low redshift bins is only 1.2cr), the 
redshift range we are probing is also not as wide 
as tha t probed by IMartini. Sivakoff. fc Mulchaevl 
(|2009h . They measure the fraction of galaxies in 
clusters that are X-ray AGN using two different 
redshift binning schemes with two and three red- 
shift bins, and restricting their sample to the two 
highest redshift bins of the three bin scheme brings 
it into better agreement with our redshift distribu- 
tion (these bins have 0.3 < z < 0.6 and z > 0.6). 
The mean X-ray AGN fraction they measure for 
these bins are 0.0031 and 0.0147, so that the in- 
crease is about a factor of 4.7. Such a large in- 
crease is unlikely in the radio source population in 
our sample. So even if qualitatively there may be 
some similar evolution (at quite low significance) 
in the radio population as the X-ray population, 
it appears to be less strong. Thus X-ray AGN 
and radio-loud AGN are likely two different AGN 
populations, and they show different redshift evo- 



lution in galaxy clusters. 

5. Conclusions 

We constrain the evolution of the bright cen- 
tral radio source population in galaxy clusters 
from redshift 0.35 to 0.95 by statistically match- 
ing FIRST radio sources with 618 galaxy clus- 
ters from a uniformly, optically selected sample 
(RCS1). The number of radio sources per clus- 
ter (0.14 ± 0.02 for clusters with 0.35 < z < 0.65 
and 0.10 ± 0.02 for clusters with 0.65 < z < 0.95) 
is consistent with previous studies using different, 
more local samples. Richer clusters have more ra- 
dio sources per cluster, also as expected from pre- 
vious studies. After normalizing by the clusters' 
mass (as derived from richness), we find that the 
number of radio sources per unit cluster mass is 
roughly constant with mass and with redshift (the 
significance of the difference between high and low 
redshift clusters is 1.5er). The number of radio 
sources (with \j\aghz > 4.1 x 10 24 W Hz -1 ) per 
10 14 M of cluster mass is 0.031±0.004. 

We also visually inspect the radio source emis- 
sion in the galaxy clusters with radio sources in 
our sample. Identifying BCGs and matching radio 
emission with optical objects, we find some (2-cr) 
evidence that the number of radio-loud BCGs is 
higher for high richness, low redshift clusters than 
for either high richness, high redshift clusters or 
low richness (at high or low redshift) clusters. If 
there is a connection between the radio activity of 
the BCG of a galaxy cluster and amount of cool- 
ing in the surrounding ICM and if the lowest red- 
shift, richest clusters are most likely to have large 
amounts of cooling gas in their inner regions, then 
we might expect to observe more radio-loud BCGs 
in nearby rich clusters. The low and intermedi- 
ate richness bins do not show any evolution in the 
prevalence of radio-loud BCGs. 

Studies of the cluster X-ray AGN population 
find strong evolution over a similar redshift range 
as we probe, while we do not observe such evolu- 
tion in the radio AGN population. This indicates 
that X-ray AGN and radio-loud AGN likely probe 
different AGN populations. 

The RCS collaboration is currently working on 
a second, larger cluster catalog based on imaging 
almost 1000 square degrees in three filters (g', r', 
z') with the CFHT MegaCam. We plan on ex- 
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tending this statistical analysis of the radio source 
population in galaxy clusters when the new sur- 
vey becomes available in order to further probe 
the effects of galaxy cluster properties on the ra- 
dio source population and to further characterize 
the cluster radio sources. 
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ing Research Council of Canada and the Canada 
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Canada-France-Hawaii Telescope (CFHT) which 
is operated by the National Research Council of 
Canada, the Institut National des Sciences de 
1'Univers of the Centre National de la Recherche 
Scientifique of France, and the University of 
Hawaii. LFB is partly supported by the Chilean 
Centro de Astrofisica FONDAP No. 15010003, 
and by CONICYT under project No. 1085286. 
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Fig. 1. — The signal-to- noise ratio of each ra- 
dial bin of the number of FIRST radio sources 
around RCS1 galaxy clusters, after a constant sur- 
face density background level has been subtracted. 
The dashed line shows the zero level at which the 
counts are equal to the background level. The ra- 
dio source positions are calculated adopting the 
clusters' redshifts for the FIRST sources. The 
background level is calculated in an annulus 5 to 
10 Mpc from the cluster center. 
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Fig. 3. — The number of radio sources per cluster, 
normalized by average cluster mass, derived from 
cluster richness, for each richness bin. The av- 
erage cluster masses are 2.1 xlO 14 M , 4.1 xlO 14 
Mq, and 1.0 xlO 15 M for the three richness bins. 
The black triangles correspond to the low redshift 
sample and the red squares correspond to the high 
redshift sample. 
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Fig. 2. — The number of radio sources per cluster 
in three different richness bins (using the B ffc # 
richness parameter). 
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Fig. 4. — Histogram of redshifts of the galaxy clus- 
ters in this study. 
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Fig. 5. — Histogram of B ffCj ^ richness parameters 
of the galaxy clusters in this study. 
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Fig. 7. — The fraction of galaxy clusters with 
radio-loud BCGs, as identified visually. The black 
triangles show the fraction for the low redshift 
sample, and the red squares show the fraction for 
the high redshift sample. 
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Fig. 6. — RCS1 z' band image of the galaxy cluster 
RCS1 J132655+302112, with radio contours from 
the FIRST image archive. The image is 0.5 x 0.5 
Mpc across at the cluster redshift of 0.37. The 
radio contours are at the 1 to 5 mJy levels, in steps 
of 1 mJy. This cluster is an example of a multiple- 
component extended radio source that looks like 
it is associated with the cluster BCG. 



Fig. 8. — The number of radio-loud non-BCG clus- 
ter members as identified visually, per cluster, nor- 
malized by the average mass (derived from rich- 
ness) for the sample of clusters in each bin. The 
black triangles show the fraction for the low red- 
shift sample, and the red squares show the fraction 
for the high redshift sample. 
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